Hard-core Yukawa model for two-dimensional charge stabilized colloids 
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The hyper-netted chain (HNC) and Percus-Yevick (PY) approximations are used to study the 
phase diagram of a simple hard-core Yukawa model of charge-stabilized colloidal particles in a two- 
dimensional system. We calculate the static structure factor and the pair distribution function over 
a wide range of parameters. Using the statics correlation functions we present an estimate for the 
liquid-solid phase diagram for the wide range of the parameters. 



the next step, we employ the widely used hypernetted- 
chain (HNC) approximation to study the liquid state 
correlation functions and freezing transition of charge- 
stabilized colloidal particles in two-dimensions interact- 
ing via the hard-core Yukawa potential (HCYP). Our 
basic aim is to test how well the HNC method mod- 
els the static properties of a colloidal system. Similar 
calculationsEl in 3D resulted in good agreement between 
the predictions of HNC theory and Monte Carlo simu- 
lations. To this end, we calculate the pair-correlation 
function and the static structure factor for a range of 
densities and temperatures. We discuss our results in 
comparison to other theoretical approaches and available 
experiments. 



I. INTRODUCTION 

Systems composed of small particles with sizes of 10's 
of nanometers dispersed in solvents are known as colloidal 
suspensions.0 With increasing interest in the properties of 
complex fluids and their various thermodynamic phases, 
colloidal systems are a subject of vast experimental and 
theoretical investigations^ 

The two-dimensional (2D) colloids are of special inter- 
est for several reasons. First, colloidal systems can be re- 
alized rather easily experimentally, as in air-water inter- 
faces or between two parallel glass plates, and advances in 
measurement techniques such as digital video microscop 



bring about a wealth of information on these systems .c 
Second, the nature of melting of 2D crystals is different 
than their 3D counterparts, making them a subject of 
intense theoretical and experimental investigations. 
Colloidal suspensions in the fluid phase may crystallize 
as their density increases or the range of interaction 
changes. Existence of a glassy phaseEj and re-entrant 
meltingS] are examples of various motivating derives 
behind some of the recent work. Possibilities of pho- 
tonic band-gap structures made out of colloidal crystals 
promise interesting applications .0 

In this work, we first use the HNC and PY approxi- 
mations for a comparison between the two approaches for 
the pure hard-core potential (PHCP) as a model for a col- 
loidal system. It is well known that the simple HNC ap- 
proximation, which omits the elementary diagrams, gives 
a good description of the large interparticle distance be- 
havior. On the other hand, the PY approximation yields 
better results for the short interparticle distances .Hill In 



II. MODEL AND THEORY 
A. Hard-Core Yukawa Model 

In this section, we introduce the model of two- 
dimensional charged colloidal particles. We consider the 
colloidal particles trapped in a surface energy well at an 
air-water interface modeled by a hard-core and a Yukawa 
interaction.i The pair potential V(r) between the col- 
loidal particles is given by 

oo for r < a , 

A exp[-ic(r- g )l for (1) 
r 1 

where A is the strength of the interaction, a is the hard- 
core radius, and k is the inverse of the Debye screen- 
ing length of the colloidal system. PHCP without the 
Yukawa tail is obtained by setting A = 0. The form of the 
interaction potential between charge-stabilized colloids is 
different from that between charged colloidal particles 
as inferred from the recent experiments of Marcus and 
Rice.Bi The role of long-range forces in colloidal systems 
has been emphasized by Noro et aZ.0 



V(r) = 



1 



B. HNC and PY formalism 



III. NUMERICAL RESULTS AND DISCUSSION 



Given the above potential to describe interactions be- 
tween charged colloidal particles, the evaluation of corre- 
lation functions reduces to a problem in classical liquid 
state theory!!. One of the basic quantities of interest is 
the pair-distribution function defined by 
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3 (|r-r'|) = (^^^(r-r i )«5(r'-r J ) 
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where p is the density of the liquid. The pair-distribution 
function gives the probability of finding a particle at r' 
when we fix another particle at r. Other relevant quan- 
tities are the pair-correlation function h(r) = g(r) — 1, 
and the static structure factor defined by 



S(q) = 1 



drh(r) exp (iq.r) 



(3) 



The pair-correlation function for a simple, isotropic 
fluid can be decomposed using the Ornstein-Zernike 
relatiorMfi 



h(r) = C(r) + p / dV C(|r - r'Qh(r') 



(4) 



where C(r) is called the direct correlation function and 
can be related to the structure factor through 



S(q) = l/[l-pC(q)) 



(5) 



A closure relation between h(r) and C(r) is needed to 
supplement the Ornstein-Zernike relation. There are 
some closure relations between h(r) and C(r). In the 
following we use the well-known the HNC and PY ap- 
proximations. In the HNC approximation, the closure 
relation takes the forr^ 1 



C(r) = exp [-/3V(r) + Y(r)] - 1 - Y(r) , 



(6) 



where V(r) is the interparticle pair potential and by con- 
struction Y(r) = h(r) — C(r). (3 = l/(fcsT) is the inverse 
temperature. On another hand, the PY closure relation 
is given by0 



C{r) = (e~ fiV ^ - l) (Y(r) + 1) 



(7) 



The self-consistent solution of Eqs. (0), (||) or ([?]) gives 
the information on the correlations in the liquid state 
within the HNC or PY schemes. 



In this section, we present the results of our numerical 
calculation of the HNC and PY equations for a system 
of colloidal particles interacting with HCYP. We have 
solved Eqs. (|), (g) or Eq. (g) for HNC and PY approx- 
imations, respectively. One can find solutions for these 
set of equations by numerical iteration taking the poten- 
tial defined in Eq. (|l|) as input. We have obtained the 
correlation functions for a wide range of values of the 
density p* = cr 2 p, screening length kct, and inverse tem- 
perature (3X. More specifically, we have calculated the 
pair-correlation function, the static structure factor, the 
collective mode energy, and direct correlation function 
which is related to the effective potential. We have then 
used the Hansen- Verlet@ criterion to estimate the phase 
line for liquid-solid transition. 



As a comparison, firstly, we compare the results be- 
tween the HNC and PY approximations. To show the 
differences, we compare the results for two extreme re- 
gions, low and high densities. In Fig. [j], we show the 
pair-distribution function at p* = 0.4 and 0.8 within the 
HNC and PY approximations for PHCP. One can see 
that the differences between two approaches increase by 
increasing density. It is evident that we have a sharp 
peak in r — a and the height of the peak increases as we 
increase the density. This means that the probability of 
finding a particle at r = a when we fix a particle in the 
origin increases with increasing density. One can observe 
the long-range oscillatory behavior in pair-distribution 
function which is a sign for liquid-solid phase transition. 
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FIG. 1. The pair-distribution function for PHCP in terms 
of ka at p = 0.4 and 0.8 within the HNC and PY approxima- 
tions. 
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FIG. 2. The static structure factor for PHCP in terms of ka 
at p — 0.4 and 0.8 within the HNC and PY approximations. 



Using Eq. (||) , one can find the static structure fac- 
tor. We depict in Fig. || the static structure factor at 
p* = 0.4 and 0.8 within the HNC and PY approxima- 
tions for PHCP. There is a peak around k = 2ir/a which 
increases with increasing density. A sharp peak in the 
structure factor at k — 2ir/a produces long-range oscilla- 
tory behavior with periodicity a in the pair-distribution 
function which is a signature for the long-range order or 
the solid phase. Furthermore, the height of this peak 
is very important in predicting of the liquid-solid phase 
transition. The Hansen- Verlet criterion simply states 
that the phase transition takes place when the value 
of the first peak in S(k) reaches a special constant. It 
is interesting that the difference between the value of 
the static structure factor within the two approximations 
near the first peak is small. This small difference allows 
the use of HNC results to be reasonable to predict the 
liquid-solid phase transition for PHCP. We used the re- 
cent results of Weber, Marx and Binder! for the density 
at which the phase transition takes place to determine 
the height of this peak in our results. This value in our 
results is about S max — 3.5, corresponding to p* — 0.89 
within the HNC approximation. 



We now report the results for HCYP within the HNC 
approximation. The static structure factor is shown in 
Fig. | at hxed /3A = 4 and na = 4 for different values 
of p* . We observe that with increasing density the peak 
structure in the static structure factor increases. Again, 
the increasing peak height of static structure factor shows 
the tendency toward the formation of solid phase. We 
also can observe that there is a shift in the position of 
the peak, which means in the high density regime the 
hard-core part of the potential is dominant compared to 
the Yukawa part and the position of the peak shifts to 
the position for the hard-core one. 



3 




We observe the same behavior by fixing na and p* 
and changing (3X or by fixing (3X and p*, and changing 
Kd . The above argument means that by changing the 
strength and range of the Yukawa part or the tempera- 
ture, one can reach the phase transition point. The re- 
sults at p* = 0.3, j3\ — 80 for different values of kg anc 
at p* = 0.4, (3\ = 4 for different values no are depicteo 
in Figs. | and |. 
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It is evident from Fig.Q that the position of the peak 
is different from the high density limit, which implies, for 
low density limit the Yukawa part has the major contri- 
bution in the liquid-solid phase transition. The position 
of the peak also has some information about periodicity 
in pair-distribution function. Generally, the position of 
the peak near the phase transition point for HCYP only 
depends on the value of the density and it is insensitive 
to other parameters. 
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FIG. 5. The static structure factor in terms of ka at 
p* — 0.4, (3\ = 4 and different values Ka. 



In Fig. ^ we show the pair-distribution function at fixed 
f3X = 4, Kd = 4 and different values of p* . It is interesting 
to note that we have a shift toward r — a as we increase 
the density, which again means that the hard-core part 
has major contribution in the high density limit. It is 
clear that the period of the oscillations also changes as 
the density increases. 
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FIG. 6. The pair distribution function in terms of rja at 
PX = 4 and no = 4 and different value of p* . In Fig | we show the var i a ti n of the first peak in terms 

of (3X at p — 0.4 and na — 4, 6, 7. It is observed that 
with increasing kct, the curves S mayi are going to tend to 
a constant, which means that hard-core part again has 
the major effect on this system in this region and we can 
thus observe the phase transition with changing density. 



It is possible to study the evolution of the first peak 
in the static structure factor when the physical quan- 
tities change, to extract some useful information about 
liquid-solid phase transition. We show in Fig. [?], the vari- 
ation of this peak in terms of j3X for a set of densities 
p* = 0.3, 0.4, 0.5, 0.6, 0.7 at kct = 4. It is clear from the 
figure that the Yukawa part is important only for low 
density and it is the source of liquid-solid phase transi- 
tion in the low density region. It is interesting that j3X 
at which the phase transition takes place also shifts to 
the right very rapidly with decreasing density, because 
the mean distance between two particles increases with 
reducing density and the effect of the Yukawa interaction 
decreases exponentially as the mean distance increases. 
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FIG. 8. The value of the first peak in the static structure 
factor in term of /3X at p = 0.4 and no = 4, 6, 7. 

From the previous figures and using the Hansen- Verlet 
criterion, we can obtain the liquid-solid phase transition 
line. We depict in Fig.KJ the phase line for different val- 
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ues of kg. It is interesting to note that all phase sepa- 
ration lines tend to the same point in the high density 
limit, which is reasonable. Another important point is 
that with increasing k<j the phase line tends to the line 
p* = 0.89. This line is equivalent to the liquid-solid phase 
transition line for pure hard-core potential. 

The phase diagram of colloidal systems in 2D is of par- 
ticular recent interest because of the subtleties in freezing 
or melting transition. Bladon and Franke]0 using an at- 
tractive aquare-well potential found a first order melting 
transition where the system goes into a liquid phase from 
a crystal phase. However, molecular dynamicsll simula- 
tions and recent experiments!! show solid-to-hexatic and 
hexatic-to-liquid phase transitions if there is an attrac- 
tive part to the colloidal potential. In an earlier work 
Lowenll did not find a hexatic phase in a Yukawa sys- 
tem. It appears that how the form of the interaction 
potential influences the possible phases in 2D colloidal 
systems is not entirely clear. Issues such as attractive or 
repulsive nature, hard or soft core, and long range part 
of the potential remain to be systematically investigated. 
In our analysis, we used the Hansen- Verlet criterion to 
determine the liquid-solid freezing transition which seem 
consistent with available results. 



investigated the effect of the Yukawa tail in the model 
potential by comparing our results with those for purely 
hard-core potential. The reliability of our correlation 
functions at small distances are tested against the PY 
approximation. Using the Hansen- Verlet criterion which 
uses the value of the first peak in the static structure fac- 
tor, we obtained the phase diagram showing the liquid- 
solid transition. 
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FIG. 9. The liquid-solid phase line transition at 
kit = 4, 6, 7. The dashed line represents the liquid-solid phase 
transition line for the PHCP. 

In summary, we have used in this paper the HNC the- 
ory to study the correlation functions of a 2D hard-core 
Yukawa fluid as a model for colloidal systems. We have 
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